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ABSTRACT 



Context. SV Psc is an asymptotic giant branch (AGB) star surrounded by an oxygen-rich dust envelope. The mm-CO line profile of 
the object's outflow shows a clear double-component structure. Because of the high angular resolution, mid-IR interferometry may 
give strong constraints on the origin of this composite profile. 

Aims. The aim of this work is to investigate the morphology of the environment around S V Psc using high-angular resolution inter- 
ferometry observations in the mid-IR with the Very Large Telescope MID-infrared Interferometric instrument (VLTI/MIDI). 
Methods. Interferometric data in the Af-band taken at different baseline lengths (ranging from 32-64 m) and position angles (73- 
142°) allow a study of the morphology of the circumstellar environment close to the star. The data are interpreted on the basis of 
2-dimensional, chromatic geometrical models using the fitting software tool GEM-FIND developed for this purpose. 
Results. The results favor two scenarios: (i) the presence of a highly inclined, optically thin, dusty disk surrounding the central star; 
(ii) the presence of an unresolved binary companion at a separation of 13.7^ g AU and a position angle of 121. S ^ ^ NE. The derived 
orbital period of the binary is 38. yr. This detection is in good agreement with hydrodynamic simulations showing that a close 
companion could be responsible for the entrainment of the gas and dust into a circumbinary structure. 

Key words. Techniques: interferometric - Stars: AGB and post-AGB - Stars: atmospheres - Stars: mass-loss - Infrared: stars 



1. Introduction 

Stars on the asymptotic giant branch (AGB) undergo phases of 
heavy stellar mass loss, which lead to the formation of circum- 
stellar envelopes. Early results obtained from interferometric 
maps of OH maser emission (e.g. Booth et al. 1981; Bowers & 
Johnston 1990) as well as CO radio line maps (Truong-Bach 
et al. 1991) were pointing to an overall spherically symmetric 
mass-loss. However, the shape of a significant fraction of 
the CO line profiles of AGB stars deviates significantly from 
that expected of a spherical envelope (e.g. Knapp et al. 1998; 
Olofsson et al. 2002; Winters et al. 2003). Such a deviation 
is evident in the case of the oxygen-rich star SV Psc, which 
shows one of the more extreme CO line shapes (Kerschbaum & 
Olofsson 1999; Winters et al. 2003). 

S V Psc is a semiregular variable star (SRb) located at a distance 
of -400 pc (Olofsson et al. 2002; Winters et al. 2003). For the 
distance determination the authors assumed a luminosity of 
the star of 4000 and 6OOOL , respectively. The mid-infrared 
dust emission features in the ISO/SWS data (Sloan et al. 2003) 
indicate an oxygen-rich, dusty environment. Kerschbaum & 
Olofsson (1999) and Olofsson et al. (2002) found that the CO 
7 = 1-0 and 7 = 2-1 lines show a clear double component in the 
outflow emission, where a very narrow feature (~ 1.5kms _1 for 
7 = 1-0, ~ 2.3 km s _1 for 7 = 2-1) is centered on a much broader 
one (-^.Skms^for 7=1-0, -lO.Skms- 1 for 7 = 2-1). This 
two-component feature was confirmed by Winters et al. (2003), 



* Based on observations made with ESO telescopes at La Silla 
Paranal Observatory under program IDs 082.D-0389 and 086.D-0069. 



but they reported slightly higher expansion velocities and 
mass-loss rates of 1.2 x lO~ 7 M yr~' (narrow component) and 
1 . 1 x 10~ 6 M yr~' (broad component). S V Psc was also detected 
in the SiO radio line emission observations (7 = 2-1, 5^1) by 
Gonzalez Delgado et al. (2003), where the two-component 
structure is visible again. 

Many scenarios have been proposed to explain this kind of line 
profile: (i) Knapp et al. (1998) interpret it as two consecutive 
winds with different age, expansion velocity and mass-loss rate. 
The narrow feature corresponds to a more recent mass-loss 
episode in this case. Winters et al. (2007) noted that this kind 
of two-wind structure may apply also to the O-rich AGB star 
EP Aqr, which shows CO line profiles similar to those of S V Psc. 
(ii) Bergman et al. (2000) found that the two-component CO 
line profile around the O-rich AGB star RV Boo can be caused 
by a low-mass disk in Keplerian rotation. They interpret the 
narrow feature as a spherically symmetric wind close to the 
star. The presence of a disk around this star was confirmed by 
Biller et al. (2005) using mid-IR adaptive optics images. Jura & 
Kahane (1999) argued that the narrow CO line profiles found 
in a number of red giant stars could be caused by long-lived 
reservoirs of orbiting molecular gas (note that these stars do 
not show double-component profiles). They interpret the disk 
as being caused by matter from previous mass-loss events 
that is captured by a binary system. This would mean that 
the disk is a direct consequence of a close binary component. 
The formation of a circumbinary disk was also predicted by 
three-dimensional hydrostatical models of Mastrodemos & 
Morris (1999). However, Knapp et al. (1998) claim that the 
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Fig. 1. iV-band spectrally dispersed Hv-coverage of the MIDI ob- 
servations of S V Psc (see Table [TJ. Contour levels range from 
8-12.5 yum (light to dark, respectively) with a step size of 0.5 /im. 

frequency of double winds appears to depend on the stellar 
chemistry and variable type and is therefore unlikely to be 
the result of binarity. (iii) Kahane & Jura (1996) suggested a 
very slow, spherically expanding wind in combination with a 
bipolar outflow for the AGB star X Her, where the latter would 
be causing the broad plateau emission. Nakashima (2005) was 
able to confirm that the broad component in the CO line profile 
could indeed be caused by a bipolar outflow and the narrow 
component is caused by a rotating disk (from interferometric 
observations of the CO J = 1-0 line). Castro-Carrizo et al. 
(2010) also found evidence for the bipolar outflow in the 12 CO 
line emission, but did not detect the Keplerian disk that was 
suggested by Nakashima (2005). 

The origin of these composite profiles, even though extensively 
studied, is still puzzling. High-angular-resolution observations 
in the mid-IR provide the possibility to study the morphology 
of the close circumstellar environment, hence providing con- 
straints on the mechanism responsible for such a line profile. In 
this paper, we present the first mid-IR spectro-interferometric 
VLTI/MIDI observations of the AGB star S V Psc, and attempt 
to find the most probable among the proposed scenarios. 

In Sect. [2] an overview of the observations of SVPsc as 
well as on the data reduction is given. The geometrical models 
and their interpretation are presented in Sect. [3] A discussion of 
the results and the conclusions are given in Sect. |4]and|5] 

2. Observations 

SV Psc was observed in 2008 (Prog. 082.D-0389) and 
2010/2011 (Prog. 086.D-0069) with the 1.8 m Auxiliary 
Telescopes of the Very Large Telescope Interferometer MIDI. 
MIDI provides spectrally resolved visibilities, photometry and 
differential phases in the iV-band (8-13yum). The baselines used 
range from 30-90 m with position angles ranging from 36-144 ° . 
The different baseline lengths provide the possibility to study 
the dusty envelope at different scales. The angular coverage (see 
Fig-[D> on the other hand, gives information on the morphology 
of the object. 



The journal of available MIDI observations is given in Table Q] 
All science and calibrator observations (with the exception of the 
calibrator YY Psc that was observed in SCI-PHOT mode) have 
been carried out in HIGH-SENS mode and in low spectral reso- 
lution (R = 30). 

In Sect. 12. II the data reduction procedure and the determination 
of errors are described. A detailed discussion of data quality and 
different selection criteria is given in Sect. 12. 21 The iV-band spec- 
trometric variability of SV Psc is discussed in Sect. 12. 31 



2.1. Data reduction 

Data were reduced using the software packages^ MIA+EWS 
1.7.1 (Jaffe 2004, Ratzka 2005, Leinert 2004). The methods used 
correspond to a coherent (EWS) and incoherent (MIA) analysis. 
The coherent analysis implements an off-line fringe tracking al- 
gorithm that aligns the interferograms before co-adding them. 
This results in a better signal-to-noise ratio of the visibility am- 
plitudes. This method also provides differential phases. In the 
incoherent method the power spectral density function of each 
scan is integrated, resulting in a squared visibility amplitude. 
The squared visibility amplitude is then integrated over time. 
The differences in the output of MIA and EWS should always 
be within 5-10% (Chesneau 2007). If this is not the case, ac- 
cording to the standard procedure, it is necessary to check if 
one of the masks used to convert the 2-dimensional frames to 1 - 
dimensional spectra is not suitable for the data reduction. EWS 
uses a fixed mask, whereas MIA creates an adaptive mask. In this 
work, data were first reduced with MIA (with MIA and EWS 
mask), and afterwards with EWS (with MIA and EWS mask), 
and finally compared with each other. 

The uniform-disk angular diameters used to derive the theoret- 
ical visibility for the calibrators as well as the corresponding 
IRAS 12/im flux and the spectral type are given in Table [2] 
The final calibrated visibilities are derived by averaging the 
MIDI visibility measurements obtained with calibrators close 
to the science target observing time. Using measurements taken 
close in time ensures that the atmospheric conditions are com- 
parable for the science and calibrator star. In this way, the trans- 
fer function should remove not only most of the instrumental 
errors, but also most of the errors introduced by the Earth's at- 
mosphere. For observations taken on 2008 Nov 03 no suitable 
calibrator was observed in HIGH-SENS mode (photometry ex- 
posures are taken after the fringe exposures) during that night. 
Therefore, one calibrator observed in SCI-PHOT mode (inter- 
ferometric and photometric signal are recorded simultaneously) 
was used to calibrate the data. As the photometry of this calibra- 
tor was also recorded independently right after the simultaneous 
recording of interferometry and photometry, we were able to re- 
duce the calibrator within the same (HIGH-SENS) mode as the 
science target. 

Error bars are derived from the standard deviation of the visibil- 
ity profiles calibrated with the given calibrators. If the derived 
error was lower than 10% or only one calibrator was available, a 
standard multiplicative error of 10% was assumed^, as the error 
introduced by the time difference in the recording of fringes and 
photometry lies between 7-15 % (Chesneau 2007). 



1 http : //www . strw . leidenuniv . nl/~ j af f e/ ews/MIA+EWS- 
Manual/index . html 

2 Except for the observations obtained on 20080 Nov 03 (see 
Sect-Hp. 
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Table 1. Journal of the MIDI Auxiliary Telescopes observations of SV Psc. 



Target 


UT date & time 


Phase [cycle] 


Config. 


Bp 


PA 


Seeing 


Mode 










[m] 


[°] 


["] 




SV Psc 


2008 Nov 03 02:17:44 


0.63 [-1] 


H0-D0 


54.5 


81 


0.44 


HIGH-SENS 


YY Psc 


2008 Nov 03 03:25:53 










0.93 


SCI-PHOT 


SV Psc 


2008 Nov 03 04:14:19 


0.63 [-1] 


H0-D0 


64.0 


73 


1.3 


HIGH-SENS 


YY Psc 


2008 Nov 03 03:25:53 










0.93 


SCI-PHOT 


SV Psc" 


2008 Oct 04 07:17:22 


0.38 [-1] 


K0-G1 


82.7 


38 


1.46 


HIGH-SENS 


SV Psc" 


2008 Oct 04 08:18:23 


0.38 [-1] 


K0-G1 


87.3 


37 


2.07 


HIGH-SENS 


SV Psc 


2008 Dec 06 01:04:13 


0.91 [-1] 


D0-G1 


65.0 


126 


0.75 


HIGH-SENS 


HD 12929 


2008 Dec 06 00:11:09 










0.76 




SV Psc 


2008 Dec 06 01:50:58 


0.91 [-1] 


D0-G1 


59.7 


125 


0.59 


HIGH-SENS 


HD 12929 


2008 Dec 06 00:11:09 










0.76 




SV Psc 


2010 Oct 09 05:34:25 


0.69 [5] 


H0-I1 


31.6 


142 


0.62 


HIGH-SENS 


aCet 


2010 Oct 09 05:14:34 










0.69 




aCet 


2010 Oct 09 05:59:17 










0.60 




SV Psc 6 


2010 Dec 08 01:57:46 


0.21 [6] 


H0-I1 


30.4 


144 


0.86 


HIGH-SENS 


SV Psc 


2011 Jan 05 01:56:59 


0.45 [6] 


I1-E0 


48.2 


92 


1.02 


HIGH-SENS 


aCet 


2011 Jan 05 01:40:14 










0.94 




aCet 


2011 Jan 05 02:14:01 










0.75 





Notes. Calibrators that were used to calibrate the data are given below the science target. The visual phase is given. The baseline configuration, 
projected baseline length and position angle of the observations are given. The observatory seeing gives the average atmospheric conditions present 
at the observatory during the observation. 
Data not used (see Sect. [2}. 



Table 2. Properties of the calibrator stars. 



HD 


Name 


Sp. T. 


Fn c 


e 










[Jy] 


[mas] 




224935 


YYPsc 


M3 III 


86.90 


7.10± 


0.07" 


18884 


o-Cet 


M1.5ina 


234.70 


11.02 


+ 0.04" 


12929 


Hamal 


K2III 


77.80 


7.43 ± 


0.52 6 



{a) http : //www . eso . org/ observing/dfo/quality/MIDI/qc/ 

calibrators_obs . html 

{b) http : //www . eso . org/ observing/etc/ 

(c) http : // simbad . u-strasbg . fr/simbad/ 
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Fig. 2. ASAS light curve of S V Psc. The grey shaded areas mark 
the range of the MIDI observations. The black square gives the 
adopted maximum for the phase determination. 



2.2. Data quality check 

During the data reduction the following selection criteria were 
applied to the science target as well as to the calibrator tar- 
gets that were used to calibrate the final visibilities, differen- 
tial phases and errors: (i) the differences in the calibrated vis- 
ibilities between MIA and EWS were within 5-10%, (ii) the 
Fourier amplitude histogram given by MIA was close to a nar- 
row Gaussian distribution, (iii) the raw visibility was smaller 
than unity, (iv) the power spectral density distribution given by 
MIA was narrow-peaked and most of the signal was located 
within the expected frequency range, and (v) the dSky parameter 
that can be adjusted by MIA (determining the region where the 
sky background is measured) had the optimal value for the given 
data set. Observations were excluded from the astrophysical in- 
terpretation if one or more of the aforementioned criteria were 
encountered. Data sets flagged with 'a' in Table[T]were excluded 
because of poor atmospheric conditions which lead to large sys- 
tematics in the calibrated measurements (Chesneau 2007). Data 
sets flagged with 'b' were not used in the interpretation because 
the differences between both reduction softwares MIA and EWS 
were larger than 10% (irrespective of the mask used). 
Observations obtained in 2008 Nov 03 show broader power 
spectral density distributions (PSD) than the others. This can in- 
troduce a systematic bias, leading to an underestimation of the 
calibrated visibility amplitude. Indeed, Fig. 4 shows that the vis- 
ibility from 2008 Nov 03 (dark grey line) is lower than that from 
2008 Dec 06 (light grey line). We therefore applied some tests 
to ensure that the change in visibility level is caused by the mor- 
phology of the star and not by the quality of the observations. 
Three pairs of observations of another science target were used. 
Each pair was observed close in time at similar baseline lengths 
and position angles with good PSD (i.e. comparable to the 2008 
Dec 06 data set) and bad PSD (i.e. comparable to the 2008 Nov 
03 data sets). The corresponding calibrated visibilities were then 
compared to each other. We found a systematic underevaluation 
of the calibrated visibilities between 5 and 25% in the case of 
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Fig. 3. Averaged MIDI flux (thin full lines) for different visual 
phases and overplotted ISO/SWS flux (thick black full line) of 
S V Psc in the A^-band. Error bars denote the standard deviation 
of all observed fluxes for the respective night and telescope. The 
ISO SWS spectrum was observed in 1998 (i.e. 33 cycles before 
the first MIDI observations). 



bad PSD. We then tested three approaches: (i) correcting for the 
visibility amplitude adding up a 25% amplitude level to the cal- 
ibrated visibilities; (ii) increasing the uncertainty level by 25% 
(see Fig. |4]i in a way that the error bars on the calibrated visibil- 
ities are ±35% (standard calibration plus systematic bias); (iii) 
not taking observations of 2008 Nov 03 into account. 
Approaches (i) and (ii) still show a difference in the visibil- 
ity amplitude between PA=73° and PA=126° (see Fig©. The 
model fitting with GEM-FIND (see Sect. 13. II ) was performed for 
all three approaches. In the case of approach (iii) the degree of 
freedom is not large enough to allow the use of the two compo- 
nent models (see Sect. l3.11 >. However, by using one component 
models we were still able to show that our observations can be 
reproduced best with elliptical models. Different two component 
as well as one component models where fitted for approaches (i) 
and (ii). The obtained results were very similar. Given this simi- 
larity, we decided to take only the second approach into account 
in the following sections. 

2.3. N-band spectrometric variability 

In Fig.|2]the light curve of SVPsc is shown (ASAS, Pojmanski 
2002). The grey shaded areas mark the range of the MIDI ob- 
servations. The visual phase of the MIDI observations is derived 
assuming a variability period of 116.33 d (taken from ACVS, 
Pojmanski 2002). The adopted phase-zero point (corresponding 
to the visual maximum of the star) for the phase determination 
is To = 2454 817 and is marked as a black square in Fig. [2] 
In Fig. [3] the MIDI spectra for different visual phases are over- 
plotted with the ISO/SWS spectrum (Sloan et al. 2003). The flux 
level of the MIDI observations is the same (within the error bars) 
as the one of ISO/SWS, implying that the level of A^-band emis- 
sion has not changed from 1998 to 201 1. 



3. Morphological Interpretation 

Figure|4] shows the MIDI calibrated visibilities for two selected 
data sets that were taken at approximately the same baseline 



Fig. 4. Spectrally dispersed MIDI visibilities of S V Psc for ob- 
servations taken at similar projected baseline lengths but differ- 
ent position angles. 



length but at different position angles. A comparison of these 
two measurements indicates a clear difference in the visibility 
level and also in the spectral shape, where the drop in the vis- 
ibility between 8 and 9.5 /vm is nearly twice as large for obser- 
vations at a position angle of 73°. As the visibility profile of a 
spherical source should be similar irrespective of the position an- 
gle, the circumstellar environment of SV Psc can be assumed to 
be not spherically symmetric. Moreover, the mid-infrared stel- 
lar spectra (see Fig. |3]l do not show any significant variation 
from phase to phase and cycle to cycle. Finally, as reported by 
Karovicova et al. (201 1), no intra-cycle/cycle-to-cycle variations 
in visibility have been detected for the oxygen-rich Mira variable 
RR Aql which presents a visual lightcurve amplitude six magni- 
tudes larger than the one of SV Psc. Therefore we can be con- 
fident that the difference in the visibility amplitude presented in 
Fig.|4]is related to a geometrical effect only. 
This difference directly proves that the spherically expanding 
wind models suggested by Knapp (1998) and Kahane & Jura 
(1996) are unlikely to be responsible for the double-component 
CO line profile in the case of SV Psc. 

In the following, the morphology of the close circumstellar envi- 
ronment of S V Psc will be constrained by fitting visibility mea- 
surements with two-dimensional, chromatic geometrical mod- 
els. 



3.1. Fitting strategy 

In order to interpret the MIDI interferometric data, the software 
GEM-FIND (GEometrical Model Fitting for INterferometric 
Data) was developed. GEM-FIND is written in the IDL lan- 
guage and allows fitting the spectrally dispersed visibility mea- 
surements with centro-symmetric and asymmetric, chromatic 
geometrical models. Each model includes a set of wavelength- 
independent and/or wavelength-dependent parameters. A de- 
tailed description of the fitting strategy and ^-minimization can 
be found in AppendixlAl 

We used different centro-symmetric and non-centro-symmetric 
models having an optically thick (i.e. one-component models) 
or optically thin (i.e. two-component models) circumstellar en- 
vironment. All spherically symmetric models show a large X^m 
implying that those models are not able to reproduce the data. 
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Fig. 5. Wavelength-dependent FWHM of the envelope's major axis and flux ratio of the central star over the envelope determined 
from the best-fitting disk model. 
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Fig. 6. Schematic view of the best-fitting models (left panel: disk model, right panel: binary model) in the sky plane at 10 /mi. 
Note that the binary component is considered to be a point source. Only one of the two possible symmetric solutions for the binary 
component location is presented here. 



Table 3. Parametric description of the best fitting geometrical 
models found in this study. 



Model 


A independent 


A dependent 


x 2 - 

rain 




fixed grid 






UD+Dirac 


Aa, AS 


J> ^prim 


0.60 


CircUD+EllGauss 


6cen ft, T] 


FWHM ma j, / 


0.49 



Notes. Aa, A6. . . angular offset of the binary component from the 
primary star; /. . . flux ratio binary/primary or central star/envelope; 
pl j m . . . the diameter of the primary component; 6 cen . . . diameter of the 
central star; ft. . . inclination angle of the ellipse; rj. . . axis ratio of the mi- 
nor/major axis of the ellipse; FWHM,„„j. . . Full Width at Half Maximum 
of the Gaussian distribution major axis; 



This is consistent with the results of Fig. [4] The optically thick 
elliptical models are also unable to provide a good fit to the 
data. This can be explained by the fact that the circumstel- 
lar environment of S V Psc is optically thin enough (i.e. M = 
1.2 x 1CT 7 M Q yr ) to allow the radiation of the central star to 
escape. An optically thin two-component model consisting of 
a central star (represented by a Uniform Disk distribution) and 
an elliptical dust envelope (represented by a Gaussian distribu- 
tion) allows a much better fit to the data. Additionally, a two- 
component model consisting of a resolved star (represented by a 
Uniform Disk distribution) and an unresolved companion (rep- 
resented by a Dirac delta function) was fitted to the data. The 
model parameters and the color-reduced minimum x 2 m i n of the 



best fitting models are given in Table[3] In the following, the two 
best fitting models are described in more detail. 

3.2. Disk model 

The circumstellar environments around low mass-loss AGB stars 
are generally considered to be optically thin (e.g. Kemper et al. 
2001). We therefore expect the stellar emission to dominate over 
the emission of the dust envelope at short mid-infrared wave- 
lenghts (i.e. ~ 8 - 9/im: see Sacuto et al. 2008). A geometrical 
model consisting of a Uniform Disk and an elliptical Gaussian 
representing a central star surrounded by a disk, respectively, 
was used. 

3.2.1. Visibility modeling 

The analytical expression of the visibility for the disk model 
can be found in Appendix lB.il This model consists of four free 
parameters: the Full-Width-at-Half-Maximum of the major axis 
FWHIvU a j, the flux ratio of the central star over the envelope fx, 
the ratio of the minor to the major axis of the ellipse r\ and the 
inclination angle of the ellipse ijj. The parameter FWHMj^ is 
chosen to be wavelength dependent to take the chromatic varia- 
tion of the opacity of the dusty environment into account. This 
results in a variation of the dimension of the structure. The pa- 
rameter f,\ is also chosen to be wavelength dependent consider- 
ing that the emission at shorter wavelengths (8-9.5 /urn) results 
from the warmer photospheric regions, whereas at longer wave- 
lengths (9.5-12.5/mi) the emission is dominated by the cooler 
dusty environment. The parameters r\ and ft are considered to 
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Fig. 7. Best-fitting disk-model visibilities (dark-grey line) superimposed on the calibrated MIDI visibilities (grey-shaded area). 
Lower right plot: Modeled visibility distribution in the wv-plane at 10pm plotted with the wv-coverage of the observations. 
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Fig. 8. Best-fitting disk-model differential phases (dark-grey line) superimposed on the calibrated MIDI differential phases (light- 
grey line). 



be wavelength independent. The diameter of the central star is 
fixed to 5.5 mas (derived from the V - K color index, van Belle 
et al., 1999) and is therefore not considered as a free parame- 
ter. The best-fitting model yields a x^ ain of 0.49 with the ma- 
jor axis of the disk inclined by an ang le of 21°+^ North-East 
and the axis ratio being O.I^qq leading to the disk-like structure 
to be seen almost edge-on. The errors on the wavelength inde- 
pendent parameters are given as the upper and lower limiting 
values of the confidence interval at a confidence level of 68%. 
The wavelength-dependent parameters are plotted in Fig. [5] The 



FWHM of the major axis of the elliptical envelope steeply in- 
creases between 8 and ~lQ/im by almost a factor of 2. Driebe 
et al. (2008) found that this signature is typical for a silicate dust 
environment. The flux ratio of the central star over the envelope, 
on the other hand, steeply decreases in this wavelength range. 
This can be explained by the fact that the emission of the central 
star becomes less significant after 9.5 fim due to the contribution 
of the silicate-rich, dusty environment. 

The left plot in Fig. [6] shows a schematic view of the structure 
for the best-fitting disk plus envelope model in the sky plane at 
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lOyum. 

The model visibilities derived with the best-fitting set of pa- 
rameters are overplotted to the observed ones in Fig. |7] The 
lower right plot in Fig. |7]shows the visibility distribution of the 
best-fitting model in the Fourier plane at lOjjm where the white 
crosses represent the Mv-coverage of the observations. 

3.2.2. Differential phase constraints 

The observed fringe phase is the visibility phase plus an 
unknown error that is introduced by atmospherical changes 
(e.g. Haniff 2007). These changes can be considered as achro- 
matic, because the refractive index of air only slowly varies with 
wavelength, i.e. errors can be related between different wave- 
lengths (Jaffe 2004, Tubbs et al. 2004). The subtraction of a 
mean phase and gradient as a function of wavenumber there- 
fore removes most of those errors. This so-called differential 
phase provides information on possible asymmetry of the source 
within the spatial and spectral region covered by the different 
interferometric measurements. For example, strong phase sig- 
nals were detected in sources harboring a circumcompanion 
disk or a circumbinary structure (Ohnaka et al. 2008, Deroo et 
al. 2007, respectively). Model phases were derived from the best- 
fitting model complex visibilities and were computed accord- 
ing to the description by Tubbs et al. (2004). They find an ac- 
curacy of the observed differential phase of 1 RMS for stable 
atmospheric conditions. However, changes in humidity during 
measurements can cause larger errors. For the measurements ob- 
tained for SV Psc the humidity and seeing levels were quite sta- 
ble throughout the night. Therefore, an error of ~5° on the mea- 
surements of the differential phases seems reasonable (Ohnaka 
et al. 2008; Deroo et al. 2007). 

Observed and modeled differential phases are plotted in Fig. [8] 
As the model is centro-symmetric, we do not expect any sig- 
nature in the differential phase. The curvature in the observed 
differential phase signatures at 54.5 m and 64.0 m is caused by 
the different indices of refraction of air present in the VLTI de- 
lay line tunnels (Deroo et al. 2007). This is normally corrected 
in the data reduction process by using the instrumental phase. 
However, in the case of the 54.5 m baseline the difference in 
time of observation of the target and the one of the calibrator 
is large, implying that the instrumental phase is not able to cor- 
rect for this curvature. On the contrary, the phase signature for 
the 48.2 m baseline seems to be real. Such a signature in the dif- 
ferential phase would point to an asymmetrical brightness distri- 
bution. A binary companion could account for this asymmetry. 

3.3. Binary model 

The binary model was constructed considering a geometrical 
model consisting of a Uniform Disk (representing the resolved 
primary component) and a Dirac delta function (representing the 
unresolved secondary component). 

3.3.1. Visibility modeling 

The analytical expression for the visibility of the binary model 
can be found in Appendix lB.2l This binary model consists of 
four parameters: the angular diameter of the primary component 
6?prim' the flux ratio of the binary component over the primary 

fx = Fbjn/Fprim and the angular coordinates of the secondary 
relative to the primary Aa and Ad. The parameter 6^ lim repre- 
sents the central star plus an extended atmosphere and is chosen 



to be wavelength dependent in order to account for the chro- 
matic variation of the opacity of the environment. This leads to a 
variation of the dimension of the AGB primary component. The 
parameter fx is also chosen to be wavelength dependent. This 
considers that the emission at shorter wavelengths (8-9.5 fim) is 
dominated by the primary star while at longer wavelengths (9.5- 
12.5 yimi) the dust material surrounding the primary star and the 
material accreted by the binary component are emitting stronger. 
The angular offset (Aa and AS) of the component, on the other 
hand, is assumed to be wavelength independent. The fitting strat- 
egy is the same as described in Appendix [A] Figure [9] shows 
the color-reduced^ 2 distribution for the geometrical model. The 
joined confidence region of the wavelength independent parame- 
ters that results from the best fit to the MIDI data is displayed by 
the white contour line. For our data a confidence level of 68% 
was used. Note that two symmetrical solutions exist, because 
the visibility determination does not allow the determination of 
the sign of the separation vector (Ratzka 2009). The best-fitting 
model yields an x-position of the binary of Aa = —29*? mas, a y- 
position of A6 = -18^j 2 mas (and for the symmetrical solution: 
Aa = 29+g mas, A<S = 18^ 2 mas) with a minimum color-reduced 
Xmin °f 0.60. Assuming a distance of 400 pc this position corre- 
sponds to a separation of 13.7^ ^ AU. The position of the bi- 
nary is well constrained by our model. A schematic view of the 
best-fitting binary model is plotted in Fig. [6] The wavelength- 
dependent parameters and fx are plotted in Fig. [10] The 
increase of the angular diameter of the primary component from 
8 to 9.5 fim can be explained by the larger opacity caused by the 
silicate dust material present in the atmosphere of O-rich AGB 
stars. The increase of the flux ratio of the binary component over 
the primary from 8 to 9.5 fim can be attributed to the fact that 
after 9.5 fim the silicate dust emission of the material that is ac- 
creted by the binary component is stronger. The visibilities ob- 
tained from the model with the best-fitting parameters are over- 
plotted on the MIDI visibilities in Fig.[TT] The lower right plot in 
Fig.Qjjshows the visibility distribution of the best-fitting model 
in the wv-plane at 10 fim. An increase of the wv-coverage using 
baselines perpendicular to the direction of the binary axis would 
help to ensure the presence of the binary companion. 

3.3.2. Differential phase constraints 

A binary system will be interpreted by an interferometer as 
an asymmetrical brightness distribution, which means that we 
expect a differential phase signature. The instrumental and at- 
mospherical corrections that are applied to the measured MIDI 
phase by the data reduction software EWS involve the removal 
of a wavenumber-dependent slope and offset from the phase 
(Jaffe 2004). The differential phase of a binary system will there- 
fore oscillate around zero (see Fig. CI in Ratzka et al. 2009 for 
an example). The missing offset and slope will prevent the ob- 
server from getting absolute astrometric information (Millour et 
al. 201 1), but will still allow to determine the position of the sec- 
ondary with respect to the primary. The differential phase is also 
the only way to disentangle between both symmetric solutions of 
the binary location. The model differential phases for both sym- 
metric solutions (dark-grey full line corresponding to the upper 
right solution in Fig. [9] and dark-grey dotted line corresponding 
to the lower left solution) are overplotted on the MIDI differ- 
ential phases in Fig. Q~2] There is no full agreement between 
the model and the differential phase measurements. However, 
because of the simplicity of the model that is used, we cannot 
exclude that a more sophisticated model, involving the full cir- 
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Fig. 10. Angular diameter of the primary and flux ratio of the binary component over the primary determined from the best fitting 
binary-model. 
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cumbinary structure, would smooth the signature revealed in the 
differential phase resulting in a better agreement with the mea- 
surements. 

The lower right plot in Fig. [T2l shows the modeled phase distri- 
bution and the observations at 10/mi in the wv-plane. According 
to this plot larger baseline lengths are needed to measure a sig- 
nificant phase signature as the central part (probed by the uv- 
coverage of the observations) is blurred by the resolved, spheri- 
cal primary component. 

4. Discussion 

The confidence region on the position of the binary component 
(see Fig. [9} leads to an angle spreading over 39.9°. The orbital 
period of the binary component was calculated by using the an- 
gular separation derived from the best fitting binary model (see 



Sect. 13 - 31 > - Assuming a mass of the primary of M pr j m =l M , a 
mass of the companion ranging from Mbi n =0.5 to 1 M Q (corre- 
sponding to the typical mass of a white dwarf or main sequence 
star) and a distance of the star of 400 pc, this leads to an orbital 
period of P or b = 38. l^ \ vr - Because our observations spread 
over 2.17 years, we expect the component to move by a max- 
imum angle of 50.4° on the plane of the sky. This makes the 
average motion of the binary companion almost consistent with 
the uncertainty we found on its position. Considering an orbital 
plane that is inclined by an angle of 32.9° relative to the plane 
of the sky, the angle of 50.4° derived from the observed orbital 
period would decrease by projection and lie within the angle of 
39.9°. 

Mastrodemos & Morris (1999) used 3D smoothed particle hy- 
drodynamics to model the influence of wind velocity, secondary 
mass and binary separation on the circumstellar environment 
of AGB stars. Considering their models 12 and 13 (Table 1 of 
the paper) that have parameters close to the ones we derived 
for SV Psc, the wind morphology would be slightly elliptical or 
bipolar. Moreover, an accretion disk around the secondary com- 
ponent of the size ranging from ~0.6 to 1 .8 AU is predicted by 
the three-dimensional models. This accretion disk would not be 
resolved with MIDI making the use of the Dirac distribution to 
model the binary component reasonable. On larger scales (few 
100 AU) gas and dust that is lost by the star is captured within the 
binary system and would form a disk-like structure. This struc- 
ture would develop from the enhancement of the density in the 
orbital plane. Such a disk could be the reason for the broad fea- 
ture visible in the CO line profile. In this case a disk in Keplerian 
rotation, as suggested by Bergmann et al. (2000) for RV Boo, 
would be a direct consequence of the presence of a binary com- 
panion. Therefore, the solution involving a disk structure seen 
edge-on does not exclude the presence of a close binary com- 
ponent. In this case MIDI would probe the most inner part (4- 
16 AU) of the rotating disk. 

In order to fully discriminate between disk-like and binary 
model, a larger wv-coverage is necessary. Additional larger base- 
lines are needed to detect a significant signal in the differential 
phase that would confirm the presence of the binary compan- 
ion. Position angles perpendicular to the binary axis would allow 
further constraining of the position of the companion, or con- 
straining the inclination angle of the disk. Moreover, additional 
measurements would allow to test more sophisticated models in- 
volving a circumbinary structure composed of the AGB star, the 
close binary component and a disk-like structure. 
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Fig. 11. Best-fitting binary-model visibilities (dark-grey line) superimposed on the calibrated MIDI visibilities (grey-shaded area). 
Lower right plot: Modeled visibility distribution in the wv-plane at 10pm plotted with the wv-coverage of the observations. 




Fig. 12. Best-fitting binary-model differential phases (dark-grey full and dotted line corresponding to upper right and lower left 
solution in Fig.|9j respectively) superimposed on the calibrated MIDI differential phases (light-grey line). Lower right plot: Modeled 
phase distribution in the Mv-plane at 10 /urn plotted with the Mv-coverage of the observations. 



5. Conclusion 

In this work we have presented the first high-angular-resolution 
mid-IR interferometric study of the close circumstellar environ- 
ment of SV Psc using VLTI/MIDI. This star is known to have 
a peculiar CO line profile that points to a non-spherical cir- 
cumstellar environment. By comparing observations taken at the 
same baseline length but at different position angles we found di- 
rect proof that the close environment of the object deviates from 
sphericity. However, the visibility level of two data sets may be 



affected by some systematic bias. Even if we derived an upper 
limit for this bias and increased the error bars on these observa- 
tions accordingly, further observations are needed to fully con- 
firm the results. 

We developed the software tool GEM-FIND allowing the fit- 
ting of centro-symmetric and asymmetric, chromatic geometri- 
cal models to the spectrally-dispersed interferometric data. We 
found two models that are able to reproduce the observed visibil- 
ities: (i) a central star surrounded by a highly inclined elliptical 
envelope having a Gaussian brightness distribution, (ii) a system 
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composed of a resolved primary star and an unresolved binary 
companion located at a separation of 14 AU. We find that the 
wavelength dependent parameters in both models are strongly 
affected by the silicate-rich, dusty environment. 
Sensitive UV and X-ray observations could give further con- 
straints on the origin and evolutionary stage of the companion. 
Futhermore, MIDI observations at larger baselines and position 
angles perpendicular to the binary axis are needed to fully dis- 
criminate between these two scenarios. A larger amount of ob- 
servations would also allow to test more complex models. For 
example, we cannot exclude the possibility that a binary com- 
panion could be the origin of a disk-like density distribution 
along the orbital plane. This has been shown by theoretical sim- 
ulations (Mastrodemos & Morris 1999) and could be tested as 
soon as additional observations are available. Such a disk could 
be the reason of the broad feature that is visible in the CO line 
profile of the star. 
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Appendix A: Detailed description of GEM-FIND 

In order to prevent solutions from falling into local minima, 
the wavelength-independent parameters included in GEM-FIND 
are varied using an equally distributed grid. This applies to all 
wavelength-independent parameters, except the diameter of the 
central star, 6* cen , which was fixed for the disk model. For each 
grid point (i.e. combination of wavelength-independent param- 
eters) a non-linear least squares fitting minimization (based on 
the Levenberg-Marquardt method: see Markwardt 2008) is per- 
formed. This allows determining the best-fitting wavelength- 
dependent parameters and a color-reduced^ 2 expressed as 



N 



for i = 1 . . .n, j = 1 ...N, 



(A.l) 



where i corresponds to the index of a given grid point (i.e. 
combination of wavelength-independent parameters) and n is the 
total number of grid points. The corresponding wavelength index 
is expressed by j, while N is the number of wavelength points. 
The number of degrees of freedom are defined as the number 
of observations minus the number of free parameters of the 
model. The best-fitting wavelength-dependent and independent 
parameters are then determined by the minimum color-reduced 
X 2 m i n - Finally, the standard uncertainties of each wavelength- 
independent parameter are derived from the minimum and max- 
imum values contained in the 68.3% confidence region denned 
by 



(A.2) 



where A depends on the number of wavelength-independent pa- 
rameters. 

Appendix B: Analytical expressions for the visibility 

B.1. Disk model 

The normalized visibilities of the two components can be ex- 
pressed as 



2Ji(n8 cen r) 



and 



exp 



'-(^FWHM^) 



2\ 



41n2 



(B.l) 



(B.2) 



The spatial frequency vector is defined as r = Vt/ 2 + v 2 . The el- 
liptical envelope is constructed by applying a rotation and com- 
pression to one axis that becomes the semi-minor axis. The cor- 
responding spatial frequency vector can be written as 



2 + ^ 



(B.3) 



where rj = a/b is the minor to major axis ratio of the elliptical 
envelope and iff is the inclination angle of the major axis counted 
from North to East. The inclined spatial frequencies can be writ- 
ten as 



- u cos iff - v sin if/ and ity = u sin iff + v cos if/ . 



(B.4) 



The final normalized visibility of the composite model can be 
written as 



V(u,v) 



|f /1 V cen ( M , v) + V* m (u, u)\ 



(B.5) 
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B.2. Binary model 

The corresponding normalized visibility of the binary model can 
be written as 

V" 1 . (u, v) + f A exp (-2ni(Aau + ASv))\ 

V(M,v) = l — — [ , (B.6) 

1 + 1 1 



where 

2Ji(n0 A - r) 

1 v nnm 7 



Vi m = ;/ nm (B.7) 



pnm 

is the normalized visibility of the primary. 
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